% %% 原始图像
% img = imread('LenaRGB.bmp');
% subplot(121)
% imshow(img);
% 
% %% 二值化
% erzhi = im2bw(img,0.4525);
% subplot(122);
% imshow(erzhi);
% 
% 
% %%
% I = imread('LenaRGB.bmp');
% subplot(1,2,1),imshow(I);
% [M,N]=size(I);
% Ic=I;
% for i = 2:M-1
%     for j=2:N-1
%         if(I(i,j)==255&I(i+1,j)==255&I(i-1,j)==255&I(i,j+1)==255&I(i,j-1)==255&I(i+1,j+1)==255&I(i+1,j-1)==255&I(i-1,j+1)==255&I(i-1,j-1)==255)
%             Ic(i,j)=0;
%         end
%     end
% end
% subplot(1,2,2),imshow(Ic);
% %%
% 
% 
% % 1 读取文件
% I=imread('LenaRGB.bmp');
% figure(1),imshow(I)
% title('原图像');
% 
% % 2 rgb图转灰度图
% if length(size(I))==3                        
%     I=rgb2gray(I);
% end     
% 
%  % 3 灰度图转二值图
% I=im2bw(I,0.4);                             
% figure(2),imshow(I)
% title('二值图');
% 
% % 4 图像膨胀处理
% B=[0 1 0;1 1 1;0 1 0];                      
% I=imdilate(I,B);
% figure(3),imshow(I);
% title('膨胀');
% 
% % 5 边缘检测（亚像素点提取）
% I=edge(I,'canny',0.9);                       
% figure(4),imshow(I)
% title('亚像素');

%%
Image = [9,2,5,4,1,1;2,17,4,16,3,3;5,3,12,4,18,3;8,20,3,9,4,6;3,3,25,5,11,8;2,3,4,6,7,8];
imshow(Image)
H2 = [-1,-1,-1;-1,8,-1;-1,-1,-1];
S = imfilter(Image,H2)
figure,imshow(S)